Understanding the Josephson current through a Kondo-correlated quantum dot 



D. J. Luitz, 1 F. F. Assaad, 1 T. Novotny, 2 C. Karrasch, 3 and V. Meden 4 'Q 

1 Institut fur Theoretische Physik und Astrophysik, Universitdt Wurzburg, 97074 Wurzburg, Germany 
2 Department of Condensed Matter Physics, Faculty of Mathematics and Physics, 
Charles University, Ke Karlovu 5, 12116 Prague, Czech Republic 
3 Department of Physics, University of California, Berkeley, California 95720, USA 

4 Institut fur Theorie der Statistischen Physik, RWTH Aachen University and 
JARA — Fundamentals of Future Information Technology, 52056 Aachen, Germany 

(Dated: January 26, 2012) 

We study the Josephson current 0-tt transition of a quantum dot tuned to the Kondo regime. The 
physics can be quantitatively captured by the numerically exact continuous time quantum Monte 
Carlo method applied to the single-impurity Anderson model with BCS superconducting leads. For 
a comparison to an experiment the tunnel couplings are determined by fitting the normal-state 
linear conductance. Excellent agreement for the dependence of the critical Josephson current on 
the level energy is achieved. For increased tunnel couplings the Kondo scale becomes comparable 
to the superconducting gap and the regime of the strongest competition between superconductivity 
and Kondo correlations is reached; we predict the gate voltage dependence of the critical current in 
this regime. 
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Introduction. Recently, hybrid superconductor-quan- 
tum dot devices have attracted much attention |T| due to 
their peculiar physical behavior determined by the inter- 
play of superconductivity of the leads and the level char- 
acteristics of the dot. Applications in nanoelectronics or 
quantum-information processing are envisaged. Among 
other properties, DC Josephson transport [2HZ] was in- 
tensively studied. Similar to the Josephson effect of or- 
dinary tunnel junctions [8 a difference <fi =/= 0, it of the 
order parameter phases of the two superconductors with 
gap A leads to an equilibrium Josephson current J run- 
ning through the system [2H7] . The focus was on carbon 
nanotube dots [5J 0HZ] with well separated single-particle 
levels, i.e., level broadening T and temperature T much 
smaller than the level spacings, simplifying the modeling 
as a single-level dot with energy e can be considered. 

It is well established both theoretically [9] as well as 
experimentally [3HZ] that the local Coulomb interaction, 
i.e. the dot charging energy U , can lead to a 0-tt transition 
of the quantum dot Josephson junction, associated to a 
first order (level-crossing) quantum phase transition from 
a singlet (0) to a doublet (tt) ground state [TO] . In fact, a 
variation of any of the system parameters U, e, A, <f>, as 
well as the tunnel couplings T L / R (with T = r^+r^j) can 
be used to tune the system across the phase boundary, if 
the others are taken from appropriate ranges. At T = 
the transition leads to a jump in J from a large and pos- 
itive (0-phase) to a small and negative value (7r-phase). 
At finite temperatures, it is smeared out and significantly 
diminished, yet the sign change of J is clearly observed 
in SQUID setups [3HH[7]. The experimental challenge in 
observing the true magnitude of the Josephson current 
to be compared with theoretical predictions consists in 
suppressing uncontrolled phase fluctuations, which can 
be achieved by using designed on-chip circuits [5] [5] . In 



such experiments J is tuned by a variation of a gate volt- 
age V g which translates into a rather controlled change 
of e [ft]. 

The physics becomes particularly interesting if the dot 
is tuned to a parameter regime in which Kondo corre- 
lations |12j become relevant for suppressed superconduc- 
tivity. It is characterized by the appearance of the Kondo 
scale (at odd dot filling) k B T K = ^TU/2exp(-nU/8T) 
|12) . Kondo physics is important if k-gT < /cbTk <C T, 
with /sb denoting the Boltzmann constant. In this regime 
perturbative methods in either U, such as self-consistent 
Hartree-Fock (HF) [l3j, or T become uncontrolled. 
Even for A ^> fceTk, at which superconductivity pre- 
vails, one expects Kondo correlations to have a signifi- 
cant impact on J. These were partly incorporated using 
a method developed for large A [13]. Other techniques 
successfully used for Kondo correlated quantum dots 
with normal leads, such as the noncrossing approxima- 
tion (NCA) [T5] , numerical renormalization group (NRG) 
[M 01], (Hirsch-Fye) quantum Monte Carlo (QMC) [T5], 
and functional renormalization group (fRG) [17] were 
extended to the present setup. With superconducting 
leads they suffer from significant conceptual or practi- 
cal limitations such as, e.g. half filling of the dot level 
(NRG) and high (NCA, QMC) or zero (fRG) temper- 
ature and, therefore, cannot be used for a quantitative 
comparison to experiments performed at temperatures 
of the order of a few tens of mK and with a wide span of 
gate voltages [FJ |5] . The regime of the strongest competi- 
tion between superconductivity and Kondo correlations 
is reached for A s» &bTk. F° r typical experimental gap 
sizes of A ps 0.1 meV [5HI]> in this regime fceTk < T is 
no longer fulfilled. Still, even for &b2k <Ta precursor 
of Kondo correlations is expected to stabilize the singlet 
phase and perturbative methods become unreliable. 
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Recently, the continuous time QMC method, called 
CT-INT in what follows, was introduced as a new tool 
to study correlated quantum dots with BCS leads [19]. 
Here, we exploit the exceptional flexibility and accuracy 
of this approach and compute J as well as the normal- 
state linear conductance G for the parameters of the ex- 
periment of Ref. [5] . Our simultaneous analysis of J and 
G reveals that the dot shows significant Kondo corre- 
lations, but superconductivity prevails as A sa 10T K . In 
the normal state it lies in the interesting and theoretically 
challenging parameter regime with k^T ss /calk ~ Unh, 
where /ish (with the Bohr magneton /zb) denotes the 
scale associated to the applied Zeeman field h used to 
destroy superconductivity. Compared to previous ap- 
proaches, we are now able to quantitatively study this 
experimentally relevant parameter regime with a numeri- 
cally exact method, and find excellent agreement between 
the experimentally measured critical current J c and the 
numerically computed one for both the 0- and 7r-phases 
(see Fig. [lj. We show that due to the fairly large left- 
right asymmetry of the tunnel couplings and the finite 
temperature the current-phase relation J{4>) is rather si- 
nusoidal even close to the 0-7r transition (see Fig. [2]), 
providing an a-posteriori justification of the extraction 
of </ c from the measured current-voltage characteristics 
of the on-chip circuits applying the extended RSJ model 
[5J[B]. Finally, using the parameters of the experiment, 
but increasing Y such that A fceTk we compute the 
gate voltage dependence of the current in the regime of 
the strongest competition between superconductivity and 
(precursors of) Kondo correlations (see Fig. p3J). 

Model and method. For the description of the single 
level quantum dot with superconducting leads we use the 
Anderson impurity model with Hamiltonian H — H dot + 
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in standard second quantized notation. In the presence 
of a Zeeman field h the single-particle energies depend 
on the orientation of the spin e a = e + g^ha [2, with 
the Lande g- factor g = 2 [2U] and a = ±1. The energy 
is shifted such that for h = 0, e = corresponds to 
the point of half-filling of the dot. The left (s = L) and 
right (s — R) superconducting leads are modeled by BCS 
Hamiltonians 
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where (without loss of generality) fa = —fa\ = <P/2. 
The quantum dot is coupled to the leads by H^ oup — 

J2k,a (*afc c Lff^ + H - c -) • Wc assume energy indepen- 
dent dot-lead hybridizations Y s = 7r^ fc |£ s fc| 2 <5(eF — tsk), 
with the Fermi energy ep. 



The CT-INT is based on an interaction expansion of 
the partition function in which all diagrams are summed 
up stochastically. The method is numerically exact 
and allows the calculation of thermodynamic observables 
with any required precision ctmc (indicated by errorbars 
in the figures) with the practical limitation that the com- 
puting time grows as l/cr^c- Details can be found in 
Ref. [19j. Here we go far beyond the proof-of-principle 
study of Ref. [TS] by considering e / 0, larger U/Y as 
well as left-right coupling asymmetries. Furthermore, we 
compute the normal-state linear conductance in the chal- 
lenging regime kp,T sa /cbTk ~ ^h. 

The Josephson current is computed as the expecta- 
tion value of the left (or right) current operator J = 

ie/hJ2ka (^Lk^ Lka d a — H.c.^. The noninteracting lead 
degrees of freedom are integrated out and one arrives 
at a formula for the Josephson current in terms of the 
imaginary-frequency Nambu-Green function Q(icu n ) of 
the dot only (directly accessible in CT-INT) [H] 



J = 2ImTr 



-Ae _I * 



G(iuj n ) 
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As our second observable we investigate the normal-state 
linear conductance G — ^ CT G a with 



G a — 
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where A a denotes the normal-state dot spectral function 
and / the Fermi function. The computation of A a from 
the (normal-state) imaginary frequency Green function 
Go(ioJn) obtained numerically by CT-INT is based on 
analytical continuation. It is found that the calculation 
of G a is much more reliable if the method detailed in 
Ref. [H] is used. As shown there the conductance can be 
written as 
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where the frequency derivative of the Green function 
has to be evaluated at imaginary frequencies iw Q which 
can differ from the Matsubara ones [22] given together 
with the weights R a in Ref. [3T]. Within CT-INT Q a 
is accessible only at the Matsubara frequencies. There- 
fore, we introduce a (real) Pade approximant Qp(oj) = 



Em - 



M-l 



Y^jLa bjUjj of degree (M, M + 1) and min- 



imize the function 



X 2 (Wi}, {h}) = ^2 {Gp(u n ) - Im[0 CT (iw„)]} 

n,m 

xCn.m {Gp(u m ) - lm[g a (iu m )}} , 

where C is the carefully bootstrapped estimate of the 
covariance matrix of the QMC data Im[C? CT (iu; n )] . The 
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degree of the Pade approximant (M, M+l) is chosen such 
that the minimal \ 2 is n °t smaller than the number of 
degrees of freedom to obtain a statistically consistent fit 
and is found to be surprisingly small with M = 3, . . . , 6. 
The Pade approximant may now be safely evaluated at 
the positions iu) a and Eq. ([5| can be used. 

Comparison to the experiment. In experiments the 
charging energy U can be determined accurately from the 
height of the Coulomb blockade diamonds obtained by 
bias spectroscopy in the normal state. The same type of 
measurement in the superconducting state reveals sharp 
features at the gap position from which A can be ex- 
tracted E]. In addition, T and h are known within 
tight bounds. The parameters which are most delicate 
to determine but strongly affect J are the level width 
T and the asymmetry T L /T R . Based on this insight, 
we proceed as follows: (i) The parameters A 7 U,T and 
h are taken directly from the experiment. Those and 
the comparison of theoretical curves for the normal-state 
conductance G(e) with the experimental ones are used 
for obtaining accurate estimates of T, T^/Fr, and the 
gate conversion factor a which relates the change of e to 
a variation of the gate voltage V g according to V g = ae 
[23] . (ii) For the complete parameter set determined this 
way, we compute the Josephson current and compare to 
the measured J c . 

We focus on the most symmetric conductance dou- 
ble peak presented in Fig. 4d) of Ref. [5]. The ex- 
perimental parameter estimates with errors of approxi- 
mately 10% are U « 3meV, A 0.1 meV, T w 75 mK, 
and h « 150 mT. In Fig. [I] a) we show our best fit of 
G(e) to the experimental result from which we extract 
r = 0.27 meV, V L /T R = 9.3, and a = 0.011 V/meV. 
At fixed U the peak separation and the peak to val- 
ley ratio are determined by T while the overall height 
is given by Tl/Tr : as can be inferred from Eq. ^ (in Q a 
only r = Tl + Tr enters) . Note that T turns out to be 
significantly smaller and /Tr significantly larger than 
the values extracted in Ref. [5] based on the assump- 
tion that the dot is in the Coulomb blockade regime. 
However, our analysis allowing for Kondo correlations 
clearly reveals that those are relevant for U/F « 11.15 
and the Kondo scale fceTk ~ 8 /ieV. It is roughly an 
order of magnitude smaller than T and of the order of 
the temperature (k^T — 6.5 ^eV) as well as the Zee- 
man energy {fi^,h = 8.7 /xeV). Therefore neither T nor h 
can be neglected when considering the normal state; the 
conductance is characterized by a split Kondo plateau 
(ridge) [23] , not to be mistaken with the Coulomb block- 
ade peaks which would be located at larger energies 
e w ±U/2 w ±1.5 meV. As an inset we show, for il- 
lustration, the normal-state spectral function at e = 
for the extracted parameters obtained from analytic con- 
tinuation of CT-INT data onto the real frequency axis 
by the maximum entropy method [25) . The appearance 
of a sharp zero energy resonance is a hallmark of Kondo 




2; 



0.0 

e [meV] 

FIG. 1. (Color online) Comparison of experimental data 5 
with the numerically exact solution of the superconducting 
Anderson model, a) Best fit of the normal-state linear conduc- 
tance with applied magnetic field used for extracting values 
of r and Tl/Fr (for details see the main text), b) Measured 
critical current vs. theoretically calculated Josephson current 
at <f> = ~ (CT-INT: symbols with line; self-consistent Hartree- 
Fock: thin lines). The arrows indicate the level positions for 
which the current phase relation is presented in Fig. [2] Inset: 
Normal-state spectral function at e = 0. 



correlations [TSJ. The splitting of the Kondo resonance 
by the Zeeman field is too small to be observable on the 
scale of the plot (but present in the data). 

In the experiments 016] J c , defined as the maximum 
of | J{4>) \ over (j) € [0, %], is extracted from current- voltage 
characteristics of the on-chip circuits using an extension 
of the standard RSJ model [2Jj]. In this analysis it is as- 
sumed that J(4>) is purely sinusoidal with its maximum 
at (j) = •§• At first glance this appears to be at odds with 
what is known theoretically for the current-phase rela- 
tion of a Josephson quantum dot in the 0-phase (half- 
sinusoidal with maximum at <j) — > n) and the transition 
region (jump from J > to J < at T = 0, smeared out 
by T > 0) [141418] . However, as it was shown already in 
Ref. [5] for an effective noninteracting model, the sizable 
left-right asymmetry and the finite temperatures of the 
experimental setups imply sinusoidal currents in the 0- 
and 7r-phase apart from very narrow ranges around the 
0-7T transitions. This conclusion is confirmed by the nu- 
merically exact CT-INT in Fig. [2] where we present J{<j>) 
for the above given parameters at the level positions in- 
dicated by the arrows in Fig. [l]b) showing \J(<j> = f , e)|. 
Apparently only for e very close to the transition the <p- 
position of the maximal current \J\ deviates observably 
from ? and yet the maximal value is still very close to 
that of \J{<f> — f )|- This gives an a-posteriori justifi- 
cation of the extraction of J c using the extended RSJ 
model and allows us to focus on the current at d> — ^ 
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FIG. 2. (Color online) Josephson current-phase relation for 
the parameters of the experiment [5] at values of e indicated 
by arrows in Fig. [l] It is rather sinusoidal even very close to 
the critical value of e (1.1 meV and 0.8 meV) and the critical 
current is thus well approximated by J {(f) = f )• 



when comparing to the gate voltage dependence of the 
critical current, as done in Fig.[l]b). Without any addi- 
tional fitting parameters we achieve excellent agreement 
in both the 0- (to the left and right of the peaks) and the 
7r-phase (central region with nearly e-independent J c ). In 
addition we show \ J{4> — § , e)| obtained for the same pa- 
rameters using the HF approach [T31 Q~7] . Whereas in the 
empty dot and doubly occupied regime |e| > 2meV this 
mean-field approximation gives decent agreement with 
the exact result (CT-INT; see also Fig. [3| it apparently 
fails in the mixed valence regime and for half dot filling 
(e « 0) in which Kondo correlations are crucial. Impor- 
tant features like the smoothing of the phase transition by 
the finite temperature and the smooth crossing through 
zero of J( | ) cannot even be obtained qualitatively. This 
emphasizes that Kondo correlations are relevant even in 
the presence of prevailing superconductivity (A s» 107k) 

Increasing Considering T — 0.4 meV and 0.5 meV 
with all the other parameters fixed at the values given 
above, we finally investigate the regime A ^b^k 
of the strongest competition between superconductiv- 
ity and (precursors of) Kondo correlations. In Fig. [3] 
\J{4>= §,e)| obtained by CT-INT is compared to HF re- 
sults. Obviously, the singlet (0) phase is stabilized by the 
correlations — an effect which HF is unable to describe. 
For the largest T K (at T = 0.5) \J(<f> = §,e)| computed 
by CT-INT only shows a precursor of the 7r-phase (the 
slight suppression close to e = 0) while HF gives a spuri- 
ous 7r-phase. It would be very interesting to measure the 
gate voltage dependence of the critical current for dots 
falling into this parameter regime, which would confirm 
the predictive power of our calculations. 

Summary. We presented a thorough study of the 
Josephson-current 0-7r transition of a quantum dot in 
the Kondo regime. A quantitative agreement to the 
measured dependence of the critical current on the gate 
voltage for a dot with Kondo correlations but prevail- 



FIG. 3. (Color online) Josephson current at 4> = f for 
the parameters of the experiment (see Fig. [T]) but with in- 
creased level broadening F = 0.4 me V and 0.5 me V (instead 
of 0.27meV) and thus increased 7k. Self-consistent Hartree- 
Fock is obviously unable to describe the strong competition 
between superconductivity and Kondo correlations in this pa- 
rameter regime and leads to a spurious 7r-phase for parameters 
for which the numerical exact solution only shows a precursor 
close to half dot filling e = 0. 



ing superconductivity was achieved. This shows that our 
minimal model is sufficient to quantitatively capture the 
relevant physics and qualifies the CT-INT as a theoreti- 
cal tool with predictive power for transport properties of 
correlated quantum dots. We further studied the regime 
of the strongest competition between superconductivity 
and Kondo correlations confirming qualitatively that the 
latter stabilize the singlet state and thus the 0-phase and 
predicting quantitatively the supercurrent, which can be 
experimentally verified. 
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